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^J • Abstract 

A microscopic model for three-cluster configurations in light nuclei is presented. It uses an 
expansion in terms of Faddeev components for which the dynamic eqations are derived. The 
model is designed to investigate binary channel processes in a compound system. Gaussian and 
oscillator bases are used to expand the wave function and to represent appropriate boundary 
conditions. We study the effect of cluster polarization on ground and resonance states of ' Be, 
and on the astrophysical S-factor of the reaction 6 Li(p, 3 He) 4 He. 
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1. Introduction 

X 

Many of the light nuclei are weakly bound. Such nuclei can change their size and shape 
considerably when interacting with other nuclei. We refer to this phenomenon as cluster 
polarization. One expects cluster polarization to play a role in reactions which involve 
light nuclei with small separation energy such as the deuteron, 6 Li, 7 Li, and so on. One 
also expects the effects to be more pronounced at small energies of the colliding nuclei, 
due to the longer interaction time intervals. 

Two different methods have been used to date to take into account polarization of 
interacting clusters. The first, introduced by Tang et al. [1|2|3|4|5| 6], considers inter- 
nal monopole excitations to describe the polarization. The second, introduced by the 
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Kiev- Antwerp collaboration [7 8 I9I10I11I12I13I14I15J , is based on collective monopole and 
quadrupole polarizations of the compound nucleus. 

In this paper we introduce a new approach in which we expand the three-cluster many- 
particle wave function into Faddeev components. This approach allows us to describe 
the proper boundary conditions for both binary and three-cluster channels. We also 
introduce two different expansion schemes: a Gaussian basis to describe bound two- 
cluster subsystems, and an oscillator basis to describe the relative motion of the third 
cluster with respect to the two-cluster subsystem. The Gaussian basis reproduces the 
intricate and complicated two-cluster bound-state behavior with a limited number of 
terms, and is thus suited to describe cluster polarization. The oscillator basis on the 
other hand allows for the proper representation of the scattering boundary condition in 
the matrix form of the Schrodinger equation. We derive a set of equations for the Faddeev 
components within the Coupled Channels Formalism. 

We apply this approach to cluster polarization in 7 Be. This system exhibits a well 
determined set of bound and resonance states, and has been thoroughly studied by many 
microscopical methods |16|17|18|19|20|21|22|23|24|25j . Moreover, two reactions which are 
connected to this nucleus, 3 He (a, 7) Be and 6 Li (p, 3 He) He, are important in astro- 
physical models |26|27j . The former reaction has been extensively investigated by different 
microscopic and semi-microscopic methods and is involved in the solar neutrino problem 
|23I20I24I21I28I29(3"U] . The latter reaction is connected to the big-bang nucleosynthesis 
and determines the abundance of light elements in the universe. It has received much 
less attention in the literature. It was investigated within a three-cluster microscopic 
model in an astrophysically relevant energy range [16j , and also in a multi-configuration 
resonating group model [18] for a wide energy range. 

We model the "' Be nucleus using a many-channel cluster wave function containing both 
two-cluster and three-cluster components. As we wish to consider both the two-cluster 
components 4 He + 3 He and 6 Li + p, we use the 4 He + d + p three-cluster configuration. 

In the reaction 6 Li (p, 3 He) He one structureless subsystem, the proton, and three 
cluster subsystems 6 Li, 4 He and 3 He are involved. They are connected with the lowest 
binary channels 4 He + 3 He and 6 Li + p which define the main properties of the bound 
and some resonance states of ' Be. 

Only 1.5 MeV is necessary to split the e Li nucleus in 4 He and d. To disintegrate the 
3 He nucleus in a deuteron and a proton the total energy of 3 He has to exceed 5.5 MeV. 
To split the 4 He nucleus into a proton and 3 H already more than 20 MeV is needed. 
This leads one to expect that cluster polarization of the 6 Li and 3 He nuclei will be 
important and should be taken into account for the low-energy states of Be, and that 
the polarization of 4 He can be neglected. 

The new method presented in this paper achieves two goals: (1) it allows to study 
the polarizability of weakly bound two-cluster systems induced by an incident cluster, 
and (2) it provides a description of the resonance structure of the resulting three-cluster 
system. 



2. The three-cluster model 

2.1. Model space and Hamiltonian 

We introduce our approach for the general case of s-shell clusters, but it can in principle 
be extended to cover clusters with an arbitrary number of nucleons. It is analogous to 
the one formulated in [31] and [32]. However, in those contributions a bi-oscillator basis 
was used to study the three-cluster interaction, while we will include both Gaussian and 
oscillator basis states. 

A Gaussian basis is a multi-parameter variational basis that can reproduce complicated 
and intricate inter-cluster wave functions with few terms, thus achieving high numerical 
precision with low computational complexity. It has been considered on different occa- 
sions in microscopic calculations, and found to be very efficient for bound states, even 
for loosely bound nuclei with proton and neutron excess [3313 4 35 36 37 38J. Its drawback 
is the non-orthogonality of the basis functions that can lead to numerical instabilities. 
We will adopt this basis to represent the (weakly) bound two-cluster subsystems in the 
three-cluster description. 

The oscillator basis is suitable for the description of bound as well as scattering bound- 
ary conditions [39 40 41 42 43 , and the corresponding matrix form of the Schrodinger 
equation [44 45 is similar to the i?-matrix theory for nuclear reactions. Due to the or- 
thogonality of the basis functions it does not suffer from numerical instabilities but it 
converges more slowly than the Gaussian basis. It was shown that an acceptable precision 
for light p-shell nuclei can be achieved with 30 to 50 oscillator functions [10113146] , In 
some model situations this number can be further reduced, even down to 3 or 5 functions, 
as was shown in |47j . The calculation of matrix elements of different operators between 
oscillator functions can be done with the technique of the Generalized Coherent States 
[48110149] . which leads to recurrence relations for the matrix elements. We will consider 
this basis to describe the scattering component in the three-cluster system. 

The wave function for s-shell clusters can be written as 

* J = A{[$! (A) $ 2 (A 2 ) $ 3 (As)] 5 [fi (xi,yi) + / 2 L (x 2 ,y 2 ) + / 3 L (x 3 ,y 3 )] )' (1) 

where $ Q (A a ) is a shell- model wave function for the internal motion of cluster a (a = 
1,2,3) and f£ (x Q ,y Q ) is a Faddeev component. The first factor describes the internal 
cluster motion and has total orbital angular momentum L = because of the s-shell 
clusters, and only its total spin quantum number is indicated. The second factor rep- 
resents the relative inter-cluster motion, and is responsible for the total orbital angular 
momentum L. We consider an LS coupling scheme so that L and S couple to the total 
angular momentum J. 

It is well know that Faddeev components are very suitable for implementing the nec- 
essary boundary conditions for binary as well as for three-cluster channels [50] . 

In the Faddeev component f£ (x Q ,y Q ), x Q is the Jacobi vector proportional to the 
distance between the (3 and 7 clusters (a, (3 and 7 form a cyclic permutation of 1, 2 and 
3), while y a is the Jacobi vector connecting the a cluster to the center of mass of the /3 
and 7 clusters: 
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For each Faddeev component we use bi-spherical harmonics 

/^(x a ,y a )= J2 fL Xa ' la ' L) (x a ,y a ){Yx a (x a )Y la (y a )} LM 
x a ,i a 



(2) 
(3) 

(4) 



which lead to the four quantum numbers X a ,l a , LM . The parity of the three-cluster 
states is then determined by the partial angular momenta: x = (— ) " ", 

The radial part of the Faddeev components /„ "' °" is obtained by using products 
of Gaussian basis functions {G\ a (x Q , b Va )} and oscillator basis functions {$n a z Q (y a , b)}, 
where 
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stands for a Gaussian function. One then immediately obtains the radial part 
fi x ^ L Hx a ,y a )= J2 cl"l a . nala G Xa (x a ,b„ a )$ na i a (y*,b) 
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Introducing the Gaussian-Oscillator bi-spherical expansion in the total wave function 
* J in H]) leads to the form 

x ^{[$ 1 (Ai)$ 2 (A 2 )$ 3 (A 3 )] 5 [GAjx Q ,^J^ /Q (y Q ^)] L } J 

a X a ,l a UocTlct 

x A{[*i (A x ) $ 2 (A 2 ) $ 3 (A 3 )] s G Aq (a: a , &„J*„„i„ (y a , 6) {**«. (x a ) yj„ (y Q )} L } J 
and will subsequently be referred to using the acronym GOB. 



The microscopic hamiltonian for a three-cluster configuration can be written as 

3 

H = f + V = Y, H£ ] + % + E V<* ( 9 ) 

a— 1 a 

i.e. a sum of three single-cluster hamiltonians H a describing the internal structure of 
each cluster, and a term responsible for the inter-cluster dynamics. The latter consists 
of the kinetic energy operator for relative motion of clusters T r and the potential energy 
of the interaction between clusters. This hamiltonian can be also expressed as a sum of 
one two-cluster hamiltonian, and terms representing the interaction of the third cluster 
with the two-cluster subsystem: 

H = H® + H& +f a + J2 %■ (10) 

The terms appearing in ([9]) and (jlpp are easily expanded in terms of the particle operators: 
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The wave function of a two-cluster subsystem 

^ = A{[$ (A / s)$ 1 (A r )} S "cj> Xa (x a )Y Xa (x Q )} J " (16) 

is expanded in Gaussian cluster functions 

^ A -=£<UkS d7) 

where 

X'l-^ =A{[<S> fi {A p )<S> 1 {A~ t )] s -G Xa {x a X)YxA* a )} Ja ■ (18) 

The bound states E„ of this subsystem [a — is the ground state, a > are 

excited or pseudo-bound states), and their corresponding eigenstates <j)" ,<7 are defined 

by \ D™'°J > and can be obtained by solving the corresponding generalized eigenvalue 

problem 

N a 

E(-,«£i 2) -4 Q) P,a)D&J = (19) 



The number of terms in (|18p . and correspondingly the number of eigenstates, is chosen 
for sufficient convergence of the ground state, and depends on the two-cluster subsystem 
labeled by a. The corresponding wave functions for the two-cluster relative motion are 

N (G) 

E D tt G ^ x ^K a ) = ^f W (20) 

Because of the Pauli principle between nucleons, it is hard to unambiguously derive a 
set of Faddeev type equations for the Faddeev three-cluster amplitudes fa(x a ,y a ) in a 
fully microscopic three-cluster description. An attempt to achieve this has recently been 
proposed in [51]. In the current paper we solve the Schrodingcr equation through the 
traditional coupled channels formalism [52153] to obtain the / Q (x a ,y Q ) amplitudes. 

The dynamic equations for the three-cluster system are easily obtained by substituting 
© in the the Schrodinger equation containing the Hamiltonian ([9]) , and in order to solve 
for the expansion coefficients C„ a x , , appropriate boundary conditions have to be 
expressed in terms of the expansion basis. This is done in the next section, exploiting 
the physical relevance of the two-cluster eigenstates discussed above. 

2.2. Boundary conditions 

In this paper we focus on the energy range between the ground state of 7 Be and the 
three-cluster threshold for 4 He+d+p disintegration. We therefore only have to consider 
binary scattering and reaction channels, and can neglect three-cluster decay. Thus only 
two-cluster asymptotics need to be included in the boundary conditions. In this case 
x a ^ Ua, i-e. one cluster is at a large distance of the other two clusters, and the latter 
will constitute a bound two-cluster subsystem. 

For large values of the Jacobi vector y a , the function fa°"°" (x a , y a ) asymptotically 
factorizes as 

f (X a ,l a ;L) {x ^ ya) „ ^y) {j , a) [S^^C-) (p a y a ) _ S COjC ^L +) (p a y a )] (21) 

for continuum states and 

/<:*«.'«,;£) (a . Q) ya) „ _^ {Xa) [s c0jC >« (-i \ Pa \ y a )] (22) 

for bound states. The entrance channel is denoted by Cq, and c a refers to the current 
channel where a stands short for all necessary (A Q , l a , . . .) quantum numbers. The mo- 
mentum p a is defined by 
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and the bound state energy E„ of the two-cluster subsystem determines the threshold 
energy of the c a channel. 

This factorization of the wave function © also occurs in the expansion coefficients 

l G v\ -n i f • T ne asymptotic region in this representation is connected to large values 

of n a (for more details see [39I40I47|10| ). and there the coefficients factorize as: 
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(24) 
(25) 

(26) 



is the classical oscillator turning point corresponding to the oscillator length b, and 
yi {Poifn a ) (respectively ij)) {Pafn a )) are the familiar radial Coulomb modified in- 
coming (respectively outgoing) wave functions, normalized to unit flux (see for instance 

152). 

The equations (|24|) and ([25]) represent the boundary condition for the expansion coef- 
ficients <^ C„ x .„ ; \ for scattering and bound states in the c a binary channel. 



2.3. The Dynamic equations 



The many-channel equations of the GOB model can be solved in three stages. 
In the first step the Schrodinger equation for all two-cluster subsystems is solved. This 
is done by diagonalizing the N v x N v matrix of the two-cluster hamiltonian 
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between the cluster Gaussian functions of (|18p . The discrete set of eigenvalues E„ cor- 
respond to bound states, or to pseudo-bound states above the threshold that are artifacts 

of the diagonalization in a finite basis. The eigenstate wave function is < D\ a ^* >. This 

step has to be repeated for every value of the partial angular momentum X a considered 
in the full calculation. 

In the second step the block matrix of the total three-cluster hamiltonian 
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is transformed to the representation of two interacting clusters using the aforementioned 
eigenfunctions. One obtains 
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This new representation exhibits the correct asymptotic behavior for large values of 
n a and ng,, in the sense that off-diagonal matrix elements coupling different channels, 
decrease to zero as n a and n& tend to infinity. 



Asymptotically the matrix has a tri-diagonal form from the kinetic energy of the rela- 
tive motion of the clusters. The diagonal matrix elements represent the interaction within 
a given channel. 

The third step in our approach consists of solving the set of equations 
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taking into account the appropriate boundary conditions to obtain either scattering or 
bound state solutions. This means that the solutions have to match the conditions (|24|) 
or (|25p respectively, beyond some matching point Ni that separates the internal and 
asymptotic parts of the wave function. Thus e.g. for scattering we look for solutions fo 
the form 
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where only the internal coefficients need to be determined. For simplicity we assume that 
Ni is identical for all channels. Inserting (|32|) into (|3Tj) then leads to 
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The solution of (|33|) then provides the explicit many-channel scattering wave function. If 
the total number of binary channels is N c , there are N c - Ni + N c - N c equations for N c ■ Ni 
expansion coefficients of the internal part of the wave function, and N c ■ N c equations for 
the determination of the 5-matrix. 



3. Results and discussion 



3.1. Parameters of the calculation 



In the current calculations we consider a Minnesota nucleon-nucleon potential (MP) 
for which we take the central part from [55] , and the spin-orbital part from [55] (data set 
IV). The exchange parameter u is fixed at u = 0.956 to reproduce the relative positions 
of the 6 Li + p and 4 He + 3 He thresholds. 

To fix the oscillator bases we use a the same oscillator radius for both the He and 
deuteron clusters. We determine it by minimizing the energy of the three-cluster threshold 



4 He + d + p, and obtain a value of b = 1.311 fm. Table Q] shows a good agreement of the 
computed threshold energies of 6 Li + p and A He + d + p compared to experiment. 

Table 1 

6 Li + p and 4 He + d + p threshold, w.r.t. the 4 He + 3 He threshold (MeV). 

Threshold MP Experiment [37] 



i. 



Li + p 4.015 

A He + d + p 5.852 



4.020 
5.493 



To fix the set of Gaussian wave functions we follow the procedure of [34158] , and 
parametrize a set of widths b v with two variational and parameters oq and q as 



a q 



1,2, 



(35) 



This has been used in [53] and [55] to obtain the ground state energy of 6 He. 
3.2. Two-cluster subsystem properties 

We first elaborate on the merits of the Gaussian basis for the two-cluster subsystems. 
To confirm its rapid convergence rate, we compare in Fig. [T] the ground state energy 
of 6 Li for both the Gaussian basis (an = 1.0 fm and q = 1.8) and the oscillator basis 
(b = 1.311 fm). We have taken the latter value considered in this paper, although it is not 
necessarily the optimal choice for this particular system. One notices that convergence is 
reached with only 4 Gaussian functions, compared to more than 20 oscillator functions. 
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Fig. 1. ! 'Li ground state energy with MP as a function of the number of Gaussian (GB) and oscillator 
(OB) cluster states. 

For the ground state of 3 He in the two-cluster d + p model, a similar situation occurs, 
and 4 Gaussian functions are sufficient with the oq = 0.9 fm, q = 1.8 parametrization. 



In Table [5] we reproduce the energies of bound and pseudo-bound states of 6 Li and 
3 He obtained with four Gaussian functions. 

Table 2 

Spectrum of bound and pseudo-bound state of G Li and ^He clusters. 

a 3 He 6 Li:L = 6 Li:L = 2 



0-5.852 -1.837 2.463 

1 1.421 3.308 4.300 

2 8.765 19.812 15.307 

3 48.774 77.898 78.670 



A standard Resonating Group Method (RGM) description takes one cluster function, 
i.e. a single (oscillator) shell-model many-particle wave function, to describe the internal 
structure of the interacting clusters. In this approximation the bound state energy of 6 Li 
is 8.800 MeV, and -3.018 MeV for 3 H e, which is way above the values obtained in Table 
[21 as was to be expected from Fig. [T] 

To obtain stable results in the three-cluster model of 7 Be for its weakly bound state, 
as well as for the elastic and inelastic scattering parameters, about 100 oscillator states 
must be considered in the calculation. This guarantees the unitarity of the calculated 
S'-matrix with high precision better than 0.1%. 

In a previous A H e + 3 He two-cluster model for 7 Be within a standard RGM approach 
[333 j O stable bound and continuous results were obtained with 30 to 50 oscillator 
states. 

By taking into account cluster polarization, the GOB model allows for more spatially 
dispersed clusters. This is confirmed by calculating the root-mean-square-radius R m of 
the two interacting clusters. For 3 He, a one cluster function approach yields R m = 
1.311 fm, while R m = 1.696 fm with four Gaussian functions. For 6 Li, these values are 
respectively R m — 1.650 fm and R m = 2.288. It is therefore natural that more oscillator 
states are necessary to properly reach the asymptotic region, because of the relatively 
large distances between the clusters compared to their sizes. 

3.3. The Spectrum of 7 Be 

In Table [3] we display the energy of the 3/2~ (bound) ground state, and the energies 
and widths of the resonance states of 7 Be obtained with the MP interaction. All the 
energies are relative to the He + 3 He threshold. There is good agreement between 

Table 3 

Ground state energy and resonance parameters (E + iT) in the GOB model of 7 Be with MP interaction 
(all in MeV and relative to the 4 i?e + 3 He threshold). 

State Theory Experiment |57l 

L = l, J n =3/2~ -1.702 -1.587 

L = 3, .r = 7/2- 2.820 + i0.130 (2.983 ± 0.05) + i (0.175 ± 0.007) 

L = 3, J w = 5/2" 5.040 + il.343 (5.143 ± 0.10) + il.20 

theory and experiment for the 3/2~ ground state. It is however slightly overbound by 

10 



0.115 MeV, and the spin-orbital splitting energy is 0.16 MeV less than the experimental 
value. The energies and widths of the lowest two resonances (7/2~ and 5/2~) are very 
close to the experimental value. The positions of the ground state and resonances of "' Be 
are displayed in Fig. [2] 
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Fig. 2. Spectrum (relative to the 4 He + 3 He threshold) of the ground and resonance states of 7 Be in 
the GOB model with MP interaction. Theory (right) and experiment (left). 

We now turn to the effect of the polarization of the two-cluster subsystems 6 Li and 
3 He on the ground state energy of ' Be. We do so comparing results with ( marked "Y") 
and without (marked " N" ) polarization of the subsystem in Table @J We suppress the 
polarization by using only a single function instead of all eigenf unctions of the corre- 
sponding two-cluster hamiltonian in the calculations. This corresponds to a rigid cluster 
throughout the calculation, whereas using the full set of eigcnfunctions allows for adapt- 
ing the size and shape of the subsystem to the presence of the third cluster. One notices 
that the polarization of 6 Li has a stronger impact than that of 3 H e. 

Table 4 

Polarization (Y: included, N: suppressed) effect on the ' Be ground state energy. 



3 He e Li 


E(MeV) 


N N 


-0.971 


Y N 


-1.413 


N Y 


-1.666 


Y Y 


-1.702 



Using the same (Y,N) approach we look at the effect of cluster polarization on the 
phase shift of 4 He + 3 He elastic scattering with total momentum J n = 7/2~ (L = 3) 
in Fig. [3l Similarly, in table [5] we look at the effect on the parameters for the 7/2 _ and 
5/2~ resonance states. 
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Fig. 3. Polarization effects on the phase shift of 4 He + 3 He elastic scattering. 

The cluster polarization substantially decreases the resonance energy and width, which 
points to an increase in the effective interaction between the clusters. The same obser- 
vation has been made in [7I8I9I10I11I13I1 4 48j for collective monopole and quadrupole 
polarizations. One again notices that the polarization of the 6 Li cluster is more important 
than that of the oi 3 He cluster. The resonance properties, obtained with 6 Li polarization 
are only marginally different from those calculated with both 6 Lz and 3 H e polarizations. 

Table 5 

Polarization effects on the resonance properties of 7 Be. 



r = 7/2" r = 5/2" 
3 He 6 Li E + ir(MeV) E + ir(MeV) 
N N 4.01 + i0.51 5.86 + i2.42 

Y N 3.69 + J0.37 5.74 + i2.24 
N Y 2.87 + i0.14 5.08 + il.38 

Y Y 2.82 + J0.13 5.04 + J1.34 



3.4. 7 Be ground state properties 

The electromagnetic observables can calculated using the explicit ground state wave 
function. In Table [75] we list root-mean-square radii (proton (R p ), neutron (R n ) and mass 
{Rm)) and the quadrupole moments (proton (Q p ), neutron (Q n ) and mass (Q m )) for the 
3/2~ ground state. As one expects, the proton radius is larger than the neutron one. The 
quadrupole moment is an indicator of the deformation of the nucleus, and a negative 
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value corresponds to a prolate deformation. The differences in these moments reflect the 
pronounced cluster structure of the Be ground state. 

Table 6 

Radii and quadrupole moments of the 7 Be 3/2~ ground state. 

r GOB MCRGMpH SCRGM[25] SVM[36] SM[59] Experiment 

R p (fm) 2.457 2.74 2.41 2.342 2.53±0.03 [60] 

R n (fm) 2.263 2.50 2.31 

R m (fm) 2.375 2.36 

Q p (e-fm 2 ) -6.245 -6.4 -6.125 -6.11 -5.153 

Q„ (e-fm 2 ) -3.739 

Q m (e-fm 2 ) -9.984 

In Table [6] we compare our GOB results with those of the multi-channel RGM of Arai 
et al (MCRGM) [16]), of the single-channel RGM of Kajino et al. (SCRGM) [25]), of the 
Stochastic Variational Method (SVM) [36] and of the Shell Model (SM) calculations [59]. 

The GOB results are close to those of SVM [36], and of the multi-channel RGM [15] . 
These results have also obtained with the MP interaction. 

The single- channel SCRGM results of [25J have been obtained with the modified 
Hasegawa-Nagata potential. The ground state quadrupole moment is comparable, whereas 
the charge radius is much larger in the SCRGM. 

The Shell Model calculation of [59] has been performed with the Bonn nucleon-nucleon 
potential. Its results for the charge radius and quadrupole moments are smaller than in 
the cluster models (MCRGM, SCRGM, SVM,GOB). This can be attributed to the fact 
that the Shell Model, involving a 10M7 or 12M7 state space, confines the inter-cluster 
distances in the ' Be nucleus. In comparison, the GOB model uses a model space of 200M7 
for the inter-cluster behavior of the dominant A He + 3 He cluster configuration. 

3.5. Comparing polarization methods 

As discussed in the introduction, one can distinguish collective polarization (quadrupole 
and monopole) for the compound nucleus as a whole (QM), RGM models with monopole 
excitations of the individual clusters (MRGM), and cluster polarization as described in 
the GOB model. It is interesting to compare these different polarization methods to 
gauge their impact on the ground state energy. 

A detailed study in MRGM is only available for ' Li [61] . so we investigate the three 
methods applied to this nucleus, as it should have deformation properties comparable to 
7 Be. 

To make the comparison consistent, we use the MP interaction parameters of [61j . 
i.e. u — 1.0174 and a spin-orbital strength of 0.821 (it is 1.0 for the other calculations 
in this paper, as suggested in [56] ). With this interaction, we have performed the QM 
calculations along the lines reported in [TS] and [T3] , a GOB calculations for a single a + 1 
channel (SGOB) and a full GOB calculation as outlined in this paper. A standard RGM 
calculation without cluster polarization (SRGM) has been included as a benchmark. 

In Table [7| the results for the different polarization models are displayed. One notices 
that the collective QM and SGOB models, which are comparable in terms of model space, 
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Table 7 

Ground state energy of 7 Li without polarization (Standard RGM) and with different types of polarization 

(see text) 

SRGMMRGM SGOB QM GOB 



-31.465 -32.027 -33.227 -33.777 -34.150 



improve with respect to the MRGM that includes only monopole cluster polarization. 
Clearly the full GOB calculation (including cluster polarization of 6 Li in the 6 Li + n 
channel) provides the best results. This was to be expected from the above 7 Be results. 

We notice that in the MRGM calculation, the monopole polarization oi 3 He is obtained 
by using four Gaussian functions resulting in -5.906 MeV of binding energy. In the SGOB 
model four Gaussian cluster functions are used to describe cluster polarization of 3 He, 
and yields a binding energy of -5.953 MeV. This confirms the above results. 

We can thus state that the cluster polarization of the channel subsystems described in 
the GOB model indeed plays an important role in seven-nucleon systems, and currently 
represents most prominent polarization type. 

3.6. Three-cluster geometry 

The expansion coefficients i C"\ .„ \ \ determine the three-cluster wave function 

W" 7 © and also the Faddeev components f a (x Q , y Q ) (a = 1, 2, 3) ^. We use the latter 
representation to analyze the wave function. 
The correlation function for the ground state is 

P a (x a ,y a ) = x 2 a y 2 a I \f a (x Q , y a ) | 2 dx a dy a (36) 

where the integration runs over the unit vectors x Q and y a . 

To interpret the the polarizability in terms of the three-cluster structure, we introduce 
the root-mean-square radii R a : 



and 



R a (s/ Q ) = \ d 9a / x 2 a \f a (x Q ,y Q )| 2 dx a /Af a (y a ) (37) 

Na(ya) = / dy a / |/ Q (x Q ,y Q )| 2 dx Q (38) 

They represent the root-mean-square radius of a two-cluster subsystem as a function of 
its distance to the third cluster (i.e. the distance of the centre-of-mass of the two-cluster 
system to the centre-of-mass of the third cluster). 

We introduce new coordinates r a and S Q whose norms (r a and S a correspond to the 
distances between the centers of mass of the clusters, that allow for a proper interpretation 
of the quantities (|3"6"]) and ([57)1 . They relate to the original Jacobi coordinates as: 



ApA, _ mr _ /^ tt (A /3 + A 7 ) s ^ 



Ap + Ay y A a + Ap + A 
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Fig. 4. Correlation function for the He + d He binary channel. 



M 4 




S<-u- P ( fm ) 



Fig. 5. Correlation function for the °Li + p binary channel 
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Fig. 2] displays the correlation function for the 4 He + 3 He channel and Fig. [5] for the 
6 Li + p channel. One notices from Fig. [4] that the distance between 4 He and 3 He is 
approximately 8 fm, and much larger than the separation of d and p. From Fig. [5] one 
observes that the binary cluster configuration 6 Li + p is surprisingly compact. Both the 
distance between He and d and between 6 Li and p are around 1 fm. These different 
geometric configurations can be related to the ground state energy relative to the thresh- 
olds of the corresponding binary channel. Indeed, the ground state is positioned at -1.702 
MeV from the 4 He + 3 He threshold and at -5.722 MeV with from the 6 Li+p threshold. 
This agrees with a very dispersed 4 He + 3 He and a compact e Li + p configuration. 
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Fig. 6. Dependence of the root-mean square radius R a of the two-cluster subsystems on the distance S a 
from the third cluster. 

In Fig. [S] we show the dependence of the root-mean-square radius of the two-cluster 
subsystems on the distance from the third cluster. One notices how 6 Li strongly adapts 
its size when the proton is at large distances (more than 15 fm) from 6 Li. When the 
e Li and the proton are near,the 6 Li + ground state gets compressed approximately 1.5 
times, and the 2 + excited state approximately two times. The 3 He nucleus, in a two- 
cluster configuration d + p, is strongly affected when the a-particlc is closer than 6 to 7 
fm. Note that without polarization all three curves in Fig. [6] would be horizontal lines, as 
the two-cluster subsystems then have a constant size. The figure illustrates the impact 
of cluster polarization of e Li and 3 He in the description of 7 Be. 

3.7. The Reaction 6 Li (p, 3 He) He 

The astrophysical S'-factor of the reaction 6 Li (p 3 He) He is obtained from the total 
cross section, which has been computed with four total angular momentum L — 0, 1, 2 
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E (MeV) 

Fig. 7. Astrophysical S-factor of the reaction 6 Li(p, 3 He) 4 He for the dominant L 
and without 3 He and e Li polarization. 



component with 



and 3. We have observed that zero angular momentum contribution dominates the low 
energy range < E < 1 MeV of the cross section. 

Fig. [7] displays the S'-factor of the 6 Li (p 3 He) He reaction calculated for L = 0. In 
contrast to the analysis of the ground state and resonance energies, the effects of cluster 
polarization arc less evident here. In Fig. [7] it is seen that 3 He polarization increases the 
5-factor, thus increasing the coupling between the channels. The 6 Li polarization only 
has a small effect on the S-factor. 

In Fig. [8] we compare the results for the S-factor with both 6 Li and 3 He cluster 
polarization with the available experimental data 62 63 64 65 66 67 68 69 70 . Notations 
and data are taken from [71] and the web site |http://pntpm.ulb. ac.be} One observes that 
our model somewhat overestimates the S'-factor in the low energy range [68 . A similar 
result was obtained in [16j . A possible reason is the lack of tensor components in the 
MP interaction. These would couple channels with different total spin and total angular 
momentum, and thus reduce the coupling between entrance and exit channels. 

4. Conclusions 

We have introduced a fully microscopic three-cluster model in which an expansion in 
terms of Faddeev components is used. The set of equations for the Faddeev amplitudes 
is derived and solved within the Coupled Reaction Channel Formalism. The dynamics of 
the three-cluster system is described by an effective two-body nucleon-nucleon potential, 
and takes into account the Pauli exchange principle correctly. 

We have used two different expansion bases: (1) a Gaussian basis suitable for the 
description for the different two-cluster subsystems, and (2) an oscillator basis to in- 
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Fig. 8. S-factor of the reaction 6 Li (p, 3 He) He. Experimental data are taken from | 62l - MA56, 
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corporate the appropriate boundary conditions for bound and continuous states in a 
many-channel compound system. Only few Gaussian states are required to describe the 
ground state of the (bound) two-cluster subsystems, such as e Li which has a high degree 
of a + d clusterization. This limits the computational effort involved in the calculations. 
As for the three-cluster results, many more oscillator states (typically 70 to 100) are 
needed to guarantee both the convergence of the ground state, and the unitarity of the 
many channel S'-matrix to sufficient precision. 

The method has been applied to investigate cluster polarizability in the ground and 
resonance states of 7 Be for which the three-cluster configuration A H e + d + p was used. 
This provides for (i) the two binary channels 3 He+ 4 He and e Li+p, which are prominent 
in the low energy region of Be, and (ii) two bound two-cluster subsystems ( 6 Li as 
4 He + d, 3 He as d + p). The latter are modeled with a Gaussian basis expansion, and 
allow to study the relative behavior of the containing clusters when e Li (respectively 
3 He) collides with a proton (respectively an a-particle). We refer to this behavior as 
"cluster polarization". The inclusion of cluster polarization in 7 Be leads to a strong 
decrease of the energy of the ground and resonance states, and reduces the resonance 
widths two to four times. The e Li polarization is more important in this respect than 
the 3 He one. 

The effect of cluster polarizability was also studied in the reaction 6 Li (p, 3 He) He by 
considering the S-factor, for which the 3 He polarization was seen to be more important. 
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